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ABSTRACT 


Resrearch conducted under this NASA/Ames-UCSC interchange has 
been concerned with numerical simulation of the structure and evo- 
lution of interstellar clouds. Steps have been taken toward an 
integrated treatment of the dynamical, thermal, and chemical 
processes entering model calculations, and a detailed study has 
been made of radiative transfer in molecular lines to allow model 
predictions to be tested against empirical data. The calculations 
have successfully reproduced and explained several observed cloud 
proerties, including abundances of complex molecular species and the 
apparent depletion of CO in dense cores. 
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I . INTRODUCTION 


The objective of work performed under this NASA/Ames-UCSC research 
interchange has been to numerically model the structure, evolution, and 
observable properties of interstellar clouds in order to better under- 
stand the processes by which interstellar material is formed into 
stars. A large variety of physical phenomena must be considered in 
such a study to represent the diversity of clouds ranging from giant 
molecular clouds with complex structure and luminous young stars to 
small globules with faint infrared stars or no evidence of star for- 
mation. The research reported here has focused on the interrelated 
roles of dynamical, thermal, and chemical processes in time-dependent 
models of collapsing clouds. In particular, we have sought to under- 
stand: 

1) The effect of the time-dependent density and temperature on 
chemical reaction rates and molecular composition. 

2) The relationship between density and temperature structure and 
the rates of heating and cooling processes. 

3) The influence of density, temperature, velocity, and chemical 
structure on molecular emission line profiles. 


I I. METHOD 


a) Chemical modeling 

The investigations of thermal, chemical, and radiative phenomena 
all utilize detailed and self-consistent numerical models for the 
dynamical evolution of gas clouds subject to forces of gravity and 
internal pressure. A hydrodynamic code provides time dependent models 
of internal structure for clouds evolving from initially uniform and 
quiescent states to configurations with large gradients of velocity, 
density and temperature. The chemical studies have used spherical 
collapse models calculated with the code of Bodenheimer 1968. (This 
one-dimensional treatment neglects effects of rotation and magnetic 
fields.) The only major modification to this code has been the in 
treatment of the gas kinetic temperature — the principal link between 
dynamics and chemistry. The temperature has been prescribed as a 
function of density and optical depth by the semi-empirical formula 
of Tarafdar, Prasad, Huntress, Villere, and Black 1985, which contains 
the essential physics of heating and cooling processes and is in 
agreement with temperatures measured in both diffuse and dense clouds. 
The density and temperature from the hydrodynamic calculations is used 
with the evolutionary chemistry code of Prasad and Huntress 1980, which 
for the present study has been extended to include 180 chemical and 
molecular species and 1626 reactions. The resulting molecular concen- 
trations and column densities are compared with cloud observations in 
Section III. 


b) Thermal structure 

A major step toward a fully consistent treatment of dynamical 
and thermal processes was undertaken with the cooperation of Dr. 

S. Tarafdar of the Tata Institute for Fundamental Research, Bombay. 
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A code was developed to follow the dynamical evolution of spherically 
symmetric clouds in which temperature is determined through heating 
by interstellar UV radiation, cosmic rays, molecule formation, and 
gravitational compression and cooling by molecular line emission and 
grain thermal emission. Because this code does not yet contain 
detailed molecular chemistry, initial applications were restricted to 
the early phases of cloud collapse. 


c) Molecular lines 

Emission line profiles of CO isotopes have been calculated for 
both rotating and nonrotating clouds modeled with the 2D hydrodynamic 
code of Black and Bodenheimer 1975. This code also has been modified 
to include the temperature parameterization of Tarafdar et al . , and 
the nonrotating models are consistent with models obtained with the ID 
code. For the dynamical calculations the gas composition is entirely 
H a , but a small admixture of CO isotopes is assumed for the emission 
profiles. Spatially uniform abundances of l2 -C0, ' 3 C0, and C isr 0 are 

used, except that at small optical depth ( / tV,‘ s in any direction £2) 
the abundance of each isotope is decreased in proportion to I'll;-, , and 
for 'Vyis ~ 4, i2 C0— =» 13 CO fractionation is modeled after results by 

Danger 1977. In some models the CO abundances have been reduced at 
high density (n(H z ) > 10* ) in accordance with the observations and 
interpretation of Wootten et al. 1978 and 1980 . 

A numerical code has been developed to solve the coupled equations 
of statistical equilibrium and raditative transfer for CO rotational 
lines. In the calculation of level populations, the radiative intensity 
is obtained from the local density, temperature, and velocity by the 
Sobolev formalism (Sobolev 1960), and collisional excitation rates are 
taken from Green and Thaddeus 1976. With the populations and source 
function thus determined, emission line profiles are obtained by 
integrating the transfer equation for selected frequencies and lines of 
sight. The line shape function for these integrations is that of Doppler 
broadening at the local kinetic temperature. This procedure is an 
improvement over the standard technique of local approximation for both 
the source function and line intensity, allowing in at least a first 
approximation for effects of large or nonmonotonic gradients in temper- 
ature, density and velocity. 


III. RESULTS 


a) Chemical modeling 

Several numerical experiments were conducted to extend the work 
of Tarafdar et al. and to address the questions of: i) large C 
abundances observed in dense clouds (N(CI)/n(C0) between 0.01 and 
0.1); ii) the non-detection of 0 x in dense clouds,- and iii) failure 
of previous chemistry studies, notably those employing hydrostatic 
models, to reproduce both the C abundance and abundances of complex 
H-C-0 molecules. 

Some investigators have suggested a C/0 abundance ratio greater 
than one to explain the large C abundance, but this leads to over- 
abundance of other carbon compounds (ie. C^H, C 3 H, C 3 H Z ). In the 
present models the abundance of C relative to CO has been enhanced in 
the cloud interior through dissociation of CO by UV radiation from 
cosmic ray excitation of the Lyman and Werner bands of H 2 . The for- 
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mation of Oi has been suppressed by reducing the rate of 0 + OH — ^ 

Ox + H on the basis of recent experimental data and by channeling 0 
into H x 0 through the assumption that H^O is the primary product of 
H 3 0 + recombinations. The non-detection of H^O in most dark clouds 
can be explained if this molecule is frozen onto grains. 

Calculations with the above assumptions have been successful in 
reproducing observed abundances of 0, 0 z, H 2 0, CO, and CaH, but the 
predicted abundances of C 3 H, and C 3 H 2 .are too small. It is found, 
however, that hydrocarbon abundances may be enhanced by artificially 
slowing the collapse rate (evolutionary time steps increased by a 
factor of 10). This suggests that non-gravitational forces (rotation, 
turbulence, and/or magnetic fields) may play an important role in 
both the dynamical and chemical evolution. Additional experiments 
have shown that increasing the temperature in high density cloud 
cores to 20-50 K (as compared with ~ 10 K obtained from the Tarafdar 
et al. formula) gives better abundances of several complex species, 
including CH 3 0H, CH 3 CN, and CH 3 C0H. This result implies that 
additional heating sources, such as newly formed stars or protostars, 
may be present although not always detected. 


b) Thermal structure 

Two projects were undertaken with the integrated hydrodynamic and 
thermodynamic code. The first calculations concentrated on compari- 
son of the temperature derived from detailed heating and cooling 
models with the semi-empirical temperature of Tarafdar et al.’ The 
code includes two numerical coefficients specifying the efficiency of 
grain photoelectron heating and the depletion of heavy elements, which 
may be varied within limits of experimental and observational uncer- 
tainty. For coefficient values within the permitted range, it is found 
that the dependence of temperature on density and optical depth is in 
good agreement with predictions of the semi-empirical formula. 

Additional calculations performed to determine the minimum mass for 
gravitational collapse as a function of cloud density showed that this 
critical mass can be an order of magnitude smaller than the Jeans mass 
for isothermal clouds. These results are qualitatively consistent 
with predictions for the stability of gas spheres with poly tropic 
temperature variation (Viala and Horedt 1974), but the critical mass 
values differ because the temperature is not polytropic when heating 
and cooling processes are treated in detail. 


c) Molecular lines 

Time sequences of models used in the calculation of CO lines 
were generated as functions of three parameters: the cloud mass, 
initial density, and initial rotational velocity. Although 
quantitative details of the evolution depend on values given to 
these parameters, the sequences show qualitative similarites that 
are reflected in general features of the spectral lines. The 
nonhomologous collapse produces a steep, centrally-peaked density 
profile. In rotating models, collapse is most rapid in the direction 
of the rotation axis, and the mass distribution becomes increasingly 
flattened. Eventually centripetal effects reverse the collapse near 
the cloud center, producing an off-center density maximum (a "ring") 
in the equatorial plane (cf. Black and Bodenheimer 1976). Both the 
collapse and rotational velocities are nonmonotonic, increasing with 
radius near the cloud center but decreasing in the outer layers. The 
gas temperature increases with radius because of decreasing density 
and optical depth. 
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An example of cloud structure and line profile calculations is 
given in Figures 1 through 4, which illustrate the evolution of a 
100 M 0 cloud from an initial state of density 100 cm -3 and angular 
velocity 5 x 10 Profiles are shown for models at which the 

maximum density in the cloud is 10 3 , 10 4 , and 10 r cm~ 3 . For each model 
the J = 1-* 0 brightness temperature of the three commonly observed CO 
isotopes is shown as a function of velocity for lines of sight at the 
indicated offsets x = R/R M o.* from the cloud center. In Figures 1 
through 3 the lines of sight are in the equatorial plane, and in Figure 
4 they are in the direction of the rotation axis. Also shown for each 
model is the variation in the equatorial plane of density, collapse 
and rotational velocities, kinetic temperature (T K ), excitation 
temperature (T x ), and CO isotope abundances. 

Several features of the line profiles give evidence of the form 
of density, temperature, and velocity distributions and of the time 
since initiation of collapse. The characteristic flattening and turn- 
over of the velocity profile in these models produces a minimum in line 
intensity because low excitation foreground material obscures warmer 
material of the same velocity. That this minimum occurs at positive 
velocity is evidence that the cloud is collapsing rather than expanding. 
Rotational motion shifts the centroid of the line profile toward higher 
velocity with increasing offset x. (At positive x the direction of 
rotation is away from the observer. ) Rotation also produces the 
differing line shapes at positive and negative x. The effects of rota- 
tion are absent in Figure 4 where the view is along the rotation axis. 

In the later models, the combination of outwardly increasing kinetic 
temperature and outwardly decreasing density produces maximum excita- 
tion temperature away from the cloud center. This leads to the off- 
center maxima in line intensity seen in Figures 2 and 3. The presence 
of two maxima in excitation temperature along a line of sight through 
the cloud center gives rise to the double emission peak in Figure 2. 
Finally, the relative intensities of the CO isotopes are seen to be 
sensitive indicators of cloud density and evolutionary state. At the 
cloud center, the line intensity of 1Z C0 decreases as increasing 
density and optical depth lower the kinetic and excitation tempera- 
tures. The intensity of the unsaturated C'^0 line increases with time 
as increasing density brings T* closer to T<. The intensity of 13 CO 
shows both effects, increasing at first and then decreasing. 

For clouds models of larger mass, the line profiles are broader 
than those of Figure 4 . The larger optical depth of these clouds 
results in lower temperature and pressure and larger collapse velocity 
and also in increased saturation of the high velocity line wings. A 
similar effect results from the assumption of larger initial cloud 
density. If the initial rotation rate is reduced, the line shapes 
approach those of the rotating cloud viewed pole-on. 

Because the comparison of theoretical and observed line profiles is 
essential both for refining physical input to the models and for 
meaningful interpretation of empirical data, effort has been made 
to present the theoretical results in forms suitable for direct com- 
parison with observations. The lightly drawn curves in Figure 3 are 
profiles that would be observed if the model cloud was placed at a 
distance of 400 pc and observed with a "state-of-the-art" telescope 
having a gaussian antenna profile of FWHM = 1' and filters of width 
0.1 km/sec. Because in this case the beam is much smaller than the 
cloud diameter (22.') and the filter width is smaller than the width of 
principal line features, the line shapes are not qualitatively altered. 
For clouds of smaller size, the effects of beam smoothing will be more 
severe . 
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Effects from the addition of a microturbulent component to the 
line broadening have also been investigated, although the results 
have only qualitative significance because the Sobolov approximation 
for the excitation temperature becomes less accurate. Line shapes are 
little changed unless the turbulent velocity is of the order of the 
macroscopic fluid velocity. Figure 5 shows the profiles of Figure 2 
with the addition of spatially uniform turbulence with characteristic 
velocity 0.5 km/s. With turbulence the profiles are wider and more 
rounded, and only the l2 -C0 lines have double peaks. 

Unfortunately observers rarely publish line profiles mapped 
across a cloud, the data most useful for comparison with numerical 
models. For this reason the present theoretical results have also 
been expressed in a form more commonly seen in the literature: 
mappings of maximum line brightness, line width, and line centroid 
velocity. Examples of these mappings are shown in Figures 6 for the 
cloud model corresponding to Figure 1. 

The theoretical line profiles were also used to assess the accuracy 
of commonly used methods by which physical parameters are deduced from 
observed profiles and simplified cloud models. Particularly note- 
worthy are results pertaining to the abundance ratio n( 13 C0)/n(H x ) . 
Several investigators (eg. Wootten et al. 1978 and 1980) find that the 
relative 13 CO abundance decreases sharply with increasing total gas (or 
H z ) density and suggest that CO freezes onto dust grains. The present 
work has shown that this apparent relationship between CO density and 
total density can result from an improper comparison of observational 
data: the lines used to derive the total gas density (H z C0 and other 
molecules of high excitation energy) are formed near the cloud center, 
whereas CO lines are formed in the outer cloud. The "observed" abundance 
variation can be reproduced using n( 13 CO) at the depth where the 13 CO 
line is formed ( T V 3 = 1) and n(H z ) at the cloud center, but if both 
densities refer to «= 1 the ratio n( C0)/n(H z .) shows little variation 
with n(H z ). (See example in Figure 7.) It is thus reasonable to suspect 
large errors in many published molecular abundances and column densities. 


IV. SUMMARY AND CONCLUSIONS 


Significant steps have been taken toward a comprehensive and inte- 
grated numerical treatment of dynamical, thermal, and chemical phenomena 
m interstellar clouds and toward the use of molecular line profiles for 
meaningful comparison of cloud models and observational data. Chemical 
evolution has been studied in the context of detailed hydrodynamical 
models, and, likewise, the thermodynamic processes that influence cloud 
temperature and reaction rates have been related to chemical and dynamical 
properties. Moreover, detailed treatment of line radiative transfer has 
shown that CO profiles are sensitive to cloud structure and evolutionary 
state, thus opening opportunities for their use in discriminating among 
alternative theoretical models and in deriving cloud parameters not 
easily inferred from observation alone. It must be realized, however, 
that m spite of these advances, further work is still needed for a 
a fully consistent and complete understanding of the physical properties 
and evolution of interstellar clouds. 
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FIGURES 


Fig. la - Emission profiles of the J = 1 — > 0 rotational transitions 
of l2 -C0, 13 CO, and-C 12f O for the evolutionary model with maximum 

density n = 10 3 cm' 3 , maximum visual extinction A v = 3, and aqe 
t = 3.4 x 10 yr. Brightness temperature is shown as a function of 
velocity for lines of sight in the equatorial plane at the indicated 
offsets x = R/R from the cloud center. 

Fig... lb - Variation with R in the equatorial plane of density, 
collapse and rotational velocities, kinetic temperature, excitation 
temperature, and CO isotope abundances for the model with maximum 
density n = 10 3 cm' 3 . 

Fig. 2a - CO emission profiles for the model with maximum density 
n = 10^ cm' 3 , maximum visual extinction A v = 14, and age 4.4 x 10 6 yr. 

Fig. 2b - Variation with R in the equatorial plane of density, 
collapse and rotational velocities, kinetic temperature, excitation 
temperature, and CO isotope abundances for the model with maximum 
density n = 10^ cm* 3 . 

Fig. 3a - CO emission profiles for the model with maximum density 
n = 10 f cm' 3 , maximum visual, extinction A v = 70, and age 4.8 x 10 yr. 

Fig. 3b - Variation with R in the equatorial plane of density, 
collapse and rotational velocities, kinetic temperature, excitation 
temperature, and CO isotope abundances for the model with maximum 
density n = 10 $ cm -3 . 

Fig. 4a - CO emission profiles for the model with maximum density 
n = 10 ^ cm' 3 , viewed in the direction perpendicular to the equatorial 
plane. 

Fig. 4b - Variation with Z along the symmetry axis of density, 
collapse and rotational velocities, kinetic temperature, excitation 
temperature, and CO isotope abundances for the model with maximum 
density n = 10 ** cm'- 3 . 

Fig. 5 - CO emission profiles for the model with maximum density 
n = 10 4 cm' 3 , with the addition of uniform turbulence of character- 
istic velocity 0 . 5 km/s . 

Fig. 6a, 6b, 6c - Contours of 13 CO brightness temperature, line 
width (FWHM) , and line velocity for the model with maximum density 
n = 10 3 cm' 3 . 

Fig. 7 - Abundance ratio n( ' 3 C0)/n(H z ) observed in a variety of 
interstellar clouds and derived from models shown in Figures 1-4. 

The filled circles are measurements by Wootten et al. 1978/ open 
circles are measurements by Snell 1979. The values labeled X are 
derived from the model profiles by the method of Wootten and Snell, 
and values labeled + are from actual model abundances at = 1 . 
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